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Abstract 

A systematic method of analysing Bethe-Salpeter equation using spectral rep- 
resentation for the relativistic bound state wave function is given. This has been 
explicitly applied in the context of perturbative QCD with string tension in the 
expansion. We show that there are only a few stable bound state mesons due to 
the small "threshold mass" (constituent mass) of quarks. The asymptotic properties 
of the bound states are analytically analysed. The spectrum is derived analytically 
and compared phenomenologically. Chiral symmetry breaking and PCAC results 
are demonstrated. We make a simple minded observation to determine the size of 
the bound states as a function of the energy of the boundstate. 
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1 Introduction 



We address relativistic bound states which are due to a causal interaction kernel. In- 
vestigation of these systems is essential to understand approximate Goldstones such as 
the physical pion. Wick-Cutkosky(WC) modeipfl |2| was one such model which was in- 
vestigated in great detail, wherein they have presented a fairly general series expansion 
technique. Here we simplify and make their formalism more transparent and infact we 
find that "deeply bound states" , those whose binding energy is comparable or more than 
the rest mass energy due to a very strong interaction kernel can be understood in a simpler 
way. An important observation has been that the most general "spectral" representation 
[3| for the bound state wave function exists and is very simple to work with. WC wave 
functions are a special class of this representation. 

In this work we will study quark-antiquark(gg) bound states in the Bethe-Salpeter(BS) 
formalism in the context of the field theoretic model (aQCD) proposed in [1]. To recapit- 
ulate the essential points of this model, string tension term(p) was explicitly incorporated 
in perturbative QCD using auxiliary fields such that ultraviolet renormalisation is assured. 
The ultraviolet (UV) behaviour remains the same as in QCD. The string tension(a) van- 
ishes asymptotically in the UV limit. In this model we will be working in the leading 
approximation and g 2 N is assumed to be small for all energies where g is the QCD gauge 
coupling constant. The infrared singular confining part of the interaction is given by the 
string tension term. Our analysis is done in Minkowski (+,-,-,-) space. 

In aQCD with our approximations we have seen that the quark propagator has no 
pole and it does not have a simple pole structure [5| unlike in WC model. Consequently the 
BS equation which involves the quark propagator has more algebraic complications. Even 
then the bound state spectral representation is still valid and this enables us to perform 
analytic calculations. Qualitatively we see that quark propagator poles are missing but 
they have "threshold masses" [S| which determines the onset of the imaginary part of the 
propagator. This is a more precise notion in our model corresponding to constituent mass 
of strong interaction phenomenology. For completeness we have presented the angular 
decomposition of the wave function in detail. For brevity we have looked at single quark 
flavour system. Our analysis can equally well handle cases of more than one flavour. 

In the BS bound state description of mesons we show that even in the presence of 
string tension there are only a few number of stable mesons and this is a consequence 
of the existence of the the threshold mass. There however are many unstable (complex 
energy) bound states and we have not made any attempt to study them systematically. 

Heavy quark bound systems under certain standard assumptions do reduce to non- 
relativistic Schrodinger theory bound systems. This is alluded to briefly as it is well 
understood in the literature. As for light mesons we derive the relationship between the 
mass of the pion and the current quark masses consistent with PCAC. 
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2 Bethe-Salpeter equation 



We address the quark-antiquark bound state problem in perturbative QCD with string 
tension. As discussed in jU Ej there are three parameters in the theory, aN, g 2 N and 
1/N of which we will treat aN as a non-perturbative parameter and the latter two per- 
turbatively. The BS equation (Fig [I]) for the quark-antiquark bound state in the 1/N 
expansion sums only the ladder graphs of V exchange' (eq.([T|)) where quark antiquark 
propagators are the non-perturbative propagators obtained by summing the rainbow 'cr 
exchange' 0. 




Figure 1: Bethe-Salpeter equation summing up ladder graphs. 



,/ -x of -\ f d4k f <rN(/)({p-k),{p + k)) g 2 N<f)((p-k),(p + k)) ^ 
t(p,p) = S(ph,J ^ ^ )YS(p) (1) 

where S(p) and S(p) are the quark propagators, a is the string tension, g is the gluon- 
fermion coupling constant and N is the number of colours. In the above we have also 
included an additional g 2 N term for the following reason. It is evident in the theory 
that the leading UV behaviour is governed by g 2 . Hence to discuss the bound state UV 
behaviour we need to consider this contribution too. With our ansatz of g 2 N small, we 
only include the leading UV behaviour (There are additional g 2 N terms interfering with 
a exchange but these are subleading in the UV regime). In this work we do not consider 
a running g or a. 

The BS amplitude is decomposed in terms of 4 x 4 Dirac matrices [UJ [7j 

<P = <Ps + 1z,<Pp + l^V + 7 M 750A + ( 2 ) 

Substituting eq.(j2J) in the BS equation , we get the following decomposition for the 
scalar, pseudoscalar, vector and pseudovector amplitudes with the propagators given by 
Sip) = t(pA(p 2 ) + B(p 2 )), S(p) = i{pA{p 2 ) + B{p 2 )) 

5 = A(AAp . p + BB) J <p s - 2(ABp„ + ABp„) J 0£ (3) 



4(AB fi + AB p)J<p s + 2iAAe a ^ s p a pp ls J <t>A» (4) 
2{AA{p ll p-p-m^ +~p l p PL ) + BByJ J , 
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7 50p = 4(AAp ■ p - 55)75 j<f> P - 2(ABp li - ABp„) 75 | ^ (5) 



7/,750^ = ±{AB p-AB t>) l5 J(f> P + 2iAAe a ^ 5 p aPpl5lh J V/1 (6) 

- 2{AA(p fl ft-p-p^ ~ BB ltl ) lh J <\>\ 

The symbol / stands for the 4-d momentum integral corresponding to the sum of the 
confining and the gluon interactions. 

t r i d 4 k gN<j)((p-k),(p + k)) g 2 N<j)((p-k),(p + k)) 
J 9 J (2tt) 41 (p - k) A (p-k) 2 ' 1 ' 

In addition the tensor components are totally determined by the above components. 

Since the momentum of the bound states has to be time like, we can go to the centre 
of mass frame of the bound state wherein the total momentum vector is given by q = 
(q = M, 0,0,0) where M is the mass of the bound state. The little group is £0(3). In 
this frame the angular momentum decomposition of the BS amplitude can be done in 
terms of 3-d solid harmonics and 0(3) scalar functions of p°,p in the following manner, 

^ = ^\p\p)y jm (8) 
for i = S,P, OV, OA i.e., Scalar ,Pseudoscalar, time component of Vector and Pseudovector 
respectively. The remaining 3 — d vector components are decomposed as [S] 

= (p^ + J^+i(pxJ)^)y jm (9) 

= z jmS <s>r (5) Y« +s ^ m (p) 



0? = (p^AP+J®2A + i(P><j)®3A)yjrn (10) 

where 5 = 0, 1, 0^ = (0 O y, 4>v) , <Pa = (4>oa, 4>a), Vjm = \p\ J Y jm , where y jm are the solid 
harmonics, J = px Vp- and Y^ +S,1 >^ m {p) are the 3-d vector spherical harmonics( 8 ). 
The relation between $|y and is given by [S] 



= SS* r<0 ' (i2) 

* 3V " Uj + i)(2j + if v -^lWTTf v (13) 

similar equations apply for the Pseudovector part. 
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3 Representation of the Wave functions 



In eq.([5 |) -eq.([ll| ) we have introduced functions of p° and p. They are Lorentz scalars as 
they are defined in the rest frame of the bound state. A convenient representation is 
required to make our analysis transparent. Consider a scalar 3-point function. In general 
this is a scalar function of momenta associated with three independent Lorentz scalar 
quantities, namely p 2 ,p 2 ,q 2 with p^ — p^ = owing to momentum conservation. Any 
scalar function associated with the 3-point function is a function of these three variables. 
There exists a spectral representation for such a function, that of Deser,et.aZ, [§]. In the 
BS wave function we are in a similar situation with one of the scalar variables namely 
q 2 fixed due to an eigenvalue condition. There are many equivalent ways of representing 
such a spectral representation. We find the most convenient one is due to [3], 



<P(p 2 ,p 2 ) = fdyTda- — (14) 

Jo Js 2 (p — yq) z — a + te 

1 

dy$((p-yq) 2 ,y) 



Note that (p — yq) 2 = (1 — y)p 2 + yp 2 — y(l — y)q 2 . 

The spectral function 4>(a, y) in general is complex and range of a can be from zero to 
infinity. For a stable bound state we know from physical considerations that it has a finite 
size and for certain range of energies of the constituents this size is not infinity. The size 
of the bound state is defined by the onset of exponential fall off in co-ordinate space. This 
is possible only if the a integration range is above some positive nonvanishing quantity 
5 2 where 5 is the inverse of the size of the bound state. In general 5 many depend on 
y. Here we will take it to be the minimum possible value in the range of integration. In 
WC model[2] the BS wavefunction can be cast into the above form where S 2 (y) is fixed 
in terms of masses of the constituents and is a series in derivatives of 8{pt — 5 2 ). This 
is also a simple case of the so called Perturbation Theory Integral representation |lUj . 
Substituting this representation for the each of the scalar functions eq.(JBJ), we can do the 
loop momentum integrals by introducing the appropriate Feynman parameter integrals 
as shown in detail in J5J. It is instructive to note the following self-reproducing property 
of the solid harmonics which follows from the defining property j2], namely, V|% m (p) = 

J d 3 kF(k 2 )y jm (k +P)= y jm (p) J d 3 kF(k 2 ) (15) 

where Fik 2 ) is a sufficiently well behaved function. 
[S][V] SECTOR: 

J §f = A(AAp-p + BB)J J $ { f - 2(ABp° + ABp°) J J (16) 
+ 2(AB + AB)f [ [ 
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J = 4(ABf + ABp°) J J $® + 2 (AAp ■ p - BB - AAp p°) J J (17) 

+ 2AA{f+p°)f J J 

J $W = 4(AB + AS) 1 1 $g - 2AA(f +p°)J J (18) 

+ 2(AAp -P-BB + 2AAf) f f 



[P][A]-SECTOR: 



| $g) = 4(AAp-p- BB) J J ¥f - 2(ABf - ABp°) J J $f A (19) 



+ 2(AB - AB)f / / $°' } 



L4 



J $S = 4(ABp° - Afip ) y y + 2{AAp -p + BB- 2AAp°p°) J J $S (20) 

+ 2AA(f+p°)f 1 J <lg 
J $g = 4(AB - AS) y | $f - 2AA(f +p°) J J $S (21) 



+ 2(A4p -P + BB + 2AAf) / / $S 



THE [V][A] Mixed Sector: for j>l 

y $g = 2(AAp. P+J Bs) y y $g+2^v y y $g (22) 

y$g = 2(AAp-p-BB)J J$®+2AAq°J J $S (23) 

= 2(AAp.p-BB)J J&£ + 2AA q yJ J*® (24) 

y$g = 2(AAp-p + BB)J j¥il + 2AAq°J J $g> (25) 
Where the symbol / stands for 

y *? ) = £dy*? ) ({p-yq) 2 ,y)) (26) 

and / / stands for 

yy^ 0) = ^ Y^C^^ 1 ^ ^^ - 2/) - - 2/)) (27) 

Q 2 N f 1 f 1 rx{p~yq) 2 

+ 



(4tt) ; 



JO JO J-oo 



for i = S,P, OV, OA, 2 V, 2A. The x J in the previous equation is replaced by for 
i = IV, 3V, 1A, 3 A and the equations are written in units of = a = 1. Also we define 



These coupled integral equations essentially become four different cases as expected 
from angular momentum algebra, namely the sum of two spin 1/2 and orbital angular 
momentum I, yields total angular momentum j as 

,7=Z®-®~ = Z®(0ffil)=Z©Z-lffiZffiZ + l (28) 

This explicit decomposition manifested in eq. (jl6|) -(|23 |) is as far as we know a new result. 



4 Asymptotic Behaviour 

First we consider the behaviour of the wave function for large space like p 2 . Here the wave 
function is real and probes the short distance behaviour. The integral equation does not 
couple the UV behaviour of the wave function to the IR or intermediate regime of the 
theory. The UV behaviour of the wavefunction is determined self-consistently by the UV 
interaction alone. The leading UV behaviour of aQCD is the same as in QCD. Using the 
asymptotic behaviour of A{p 2 ) and B{p 2 ), 

A{p 2 ) ~ 1 (29) 



—p 2 y2a ln(— p 2 ) 
B(P>) ~ !2fc£ (30) 

— P A 

we adopt the same procedure as shown in 

The leading order asymptotic behaviour of the BS amplitudes in the [P] [A] and [S] [V] 
sector are the same. Its in the next to leading order(NLO) that they differ. The leading 
order BS amplitudes go like, 



For j > 1 



^ ~ (j>f 



'oa ~ <Pqv 



hA ~ 01V 



hv 03A 



(-P 2 ] 


H+ 2 (ln( 


-p 2 )) 1+ lTT 




1 




(-P 2 : 


)i+ 2 (ln( 






1 




(-p 2 . 


)J+ 3 (ln( 


-P 2 )) 1_ ^2 




1 




{-p 2 \ 


)i+ 2 (ln( 






1 





[-p 2 y+ 3 (\n(-p 2 )) 1+ m 



(31) 
(32) 
(33) 

(34) 
(35) 
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We have not used running a in the above analysis. It is seen that A(p 2 ) function 
dominates as expected and the asymptotic behaviour of the wave function is independent 
of a due to the dependence of A(p 2 ) on a as given in eq.(|'29p. It is also evident that no 
further infinite renormalisations are needed as the BS wavefunction asymptotic behaviour 
is sufficiently small that all momentum integrals are finite. This demonstrates that the 
theory is made finite by the standard wavefunction, mass and string tension renormalisa- 
tions alone. The leading aysmptotic behaviour as shown in eq.([31| ) -([35| ) is the same both 
in QCD and oQCD. For completeness we mention that if we ignore the g 2 N term in 
eq.((T]) then the asymptotic analysis, yields similar results as eq.(|3T |) -(|33 j) with the powers 
of p 2 decreased by one and all ln(— p 2 ) powers are zero. 

5 Spectrum of Light Mesons 

The BS equation is simplified considerably in their algebraic complexity. Generically it is 
very much like in the WC model. Major difference being that the propagator functions 
are complicated functions unlike simple poles in the WC model. The eigenvalue problem 
is well defined once the explicit A{p 2 ) and B{p 2 ) functions of the quark propagator are 
given. 

The most important properties of A(p 2 ) and B(p 2 ) that we exploit is that they are 
analytic functions near p 2 = and the onset of non-analyticity is near the threshold mass 
fh, i.e., p 2 = fh 2 in units of a. Considering the BS wavefunctions ()16)) - (j25|) we first note 
that these functions are real and analytic for p 2 < fh 2 and p 2 < rh 2 . This is equivalent 
to saying that for all eigenvalues q such that qo < fh + rh the wave functions are real and 
analytic. It is also evident from the standard arguments pi 0] that for qo > fh + rh the 
wavefunctions are necessarily complex and perhaps even unstable, i.e., the eigenvalue q 
itself may be complex. 

For light quarks where the renormalised mass(current mass) m is much smaller than 
a, we have shown [S] that threshold mass(constituent mas) rh 2 <C o, indeed we estimated 
that rh 2 ~ .02a. For stable mesons qo < rh + fh, consequently go y/&. Therefore all 
stable bound states in this system are necessarily deeply bound. For simplicity we ignore 
the gluon coupling and keep only the string tension contribution to the BS equation. We 
solve the BS equation at p 2 = p 2 = where we know explicitly the propagator functions. 
Consider the case when the BS wavefunction is non vanishing at p 2 = p 2 = 0, since q 2 ^ 1, 
we negelect the xq 2 dependence in the r.h.s of eq. (jl6|) - PBT) . Then all x integrations can be 
done explicitly (The y integration can be done formally on both sides). Consequently the 
BS equation reduces to an ordinary matrix eigenvalue equation in each of the different 
sectors at p 2 = p 2 = 0. Solution to these homogenous equations is guaranteed if the 
corresponding determinant of the matrix is zero. Noting that all our calculations are 
valid only if q 2 > the relevant solutions resulting from the vanishing of the determinant 
are given below. The explicit answers are given for renormalised quark mass, m = rh. 



[S][V] SECTOR: q 2 



2a 2 L(l - 2b 2 H) 



(36) 
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For small quark masses m we have for j = and 1, M 2 (j p ) where P is the intrinsic parity. 

M 2 (0+) « J + ^rn (37) 

M 2 (l-) « JL + ?lZI m (38) 
V ; 138 6348 V ; 



, (2b 2 L-l)(4b 2 L + l) 



[P][A] SECTOR : g ^ ^L(i ; 2 TO ^ 



For small quark masses m we have for j = 



M 2 = M 2 (0~) « -m (40) 



[V][A] SECTOR 



For small quark masses m we have for j 



M 2 (r 



M 2 (\-\ 



(2blH+l)(2blL-l) 
(26 ( 2 ,f/-l)(2fegL+l) 



23 121 

1 m 

160 160 

7 13 

1 m 

160 32 



(41) 
(42) 



where 



m - Vm 2 + 16 b 

bo = 5 ; a = — (43) 

8 m — o 

iJ= / -L = / rfx^ (44) 

Jo 1 — x Jo 1 — X 



Although the angular momentum j can become arbitrarily large, we find for larger j than 
what we have considered, q 2 becomes negative, thus negating our initial assumption that 
go is time like. Hence these are discarded. 

In addition we can have solutions where the BS wavefunctions vanish at p 2 = p 2 = 0. 
Since they have to be analytic they have to vanish as integer powers of p 2 or p 2 . Noting 
that in our representation eq. ()14|) . this can only be of the form {{p — yq) 2 ) n and in the limit 
of small q 2 this also approximately vanishes. Implementing these kind of wavefunctions 
we get precisely the radial excitations. Trivial algebra shows that approximately the 
eigenvalues M 2 are functions of j ' + n only. All these eigenfunctions can be given Taylor 
expansions in p 2 and p 2 just as we did in [5 for the quark propagator. This is a double 
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series and convergence properties of this series is technically more cumbersome to handle 
and we have postponed it for later study. 

Indeed the above conditions are only necessary conditions for the existence of the 
bound state. Sufficient conditions have yet to be stated. In addition to the above there 
are many complex solutions with real part of qo > m + rh. While these are acceptable as 
eigenvalue conditions, these should be truly taken as unstable resonances. 

Many of the eigenvalues both real and complex have to satisfy the finite size criterion. 
Namely the size of a bound state or the extent to which the wave function is spread out 
should be finite i.e., 5 2 in eq.()14jl should be non-zero. We are unable to estimate this 
analytically from the BS equation but a heuristic argument to be stated later suggests 
that all eigenvalues with q° < rn + rh, where rh and rh are the threshold masses, can exist. 

Let us look at the phenomenological implications of our spectrum. Figj2j gives eigen- 
value M 2 versus the quark mass m, both of which are in units of a and \f& respectively. 
We first compare pseudoscalar 0~ with the rest. This is very well understood in QCD 
From first principles both in continuum and lattice under wide circumstances one 
can show that the lightest meson is the pseudoscalar, in particular M(0~) < M(l~). This 
is also borne out by experimental data. In our theory this inequality cannot be formally 
shown to be valid, however it is maintained for renormalised mass m less than .07 and it is 
disobeyed for larger m. So we have to conclude from this that this theory is qualitatively 
different from QCD for m > 0.07. 

Next we will take M v and M p mesons as given to fix physical values for a and the mass 
of the up quark, m u (for down quark we take m d = m u ). We find from eq.pSjl and eq.(|4"U)l. 
o = (AGeV) 2 and m u = 6MeV 12. The Scalar + turns out to be very heavy(1.56GeV r ) 
and cannot certainly be a stable bound state as it is greater than 2rh u , where rn u is the 
threshold mass of the up quark which we estimate in this model to be about .56GeV. 
This implies all up and down quark bound states less than 1.2GeV are stable. 

In the strange quark sector, the 1~ vector bound state of ss is unambigously known 
to be 0(1020). From Figj2j, it can be inferred that the strange quark mass(m s ) is about 
0.26GeV and a pseudoscalar 0~ will have a mass 0.95GeV which corresponds to the r\ 
meson. That is, in this scenario most of the mass of the rj meson can be thought of 
as coming from ss (flavour mixing is not attempted in our calculations). In our model 
strange quark mass is very close to the crossover regime where 1~ and 0~ cross, beyond 
which the model qualitatively fails to be QCD like. 

Eq (|40j) demonstrates the well known consequence of PCAC namely, the square of the 
mass of the pion is proportional to the current quark mass m and the proportionality 
constant in this model is §v^ — 3Gev. Furthermore we see that there are two states way 
above the threshold, namely + and 1 + . + state is the well known "a particle". It is 
clearly extraordinarily massive and is expected to be unstable even in this lowest order 
calculation as it exceeds the threshold energy. 

Physical spectrum is expected to show a lot of mixing in flavour neutral particles. This 
can be anticipated in this model purely because the threshold mass for all the flavours is 
about the same [5], hence corrections can become dominant due to kinematical reasons 
alone . By carrying out calculation we can have better fit to phenomenology. 
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0.25 




Figure 2: Meson Spectrum with j p . The one in the bracket is a radial excitation. Only 
mesons with q 2 < 4m 2 are stable 

Now we make a semi-analytic discussion as to the size of the bound states from general 
considerations of BS equation. These equations have a generic form as 



<P = SS j K(j) (45) 



where S(p) and S(p) are the propagator functions of the constituents and K is the interac- 
tion kernel. In general spatial length scales can be present from the interaction kernel. For 
our discussion we shall assume that the interaction kernel has long range like the massless 
gluon. Even the string tension term is long ranged. For these type of kernels, the length 
scales come from the propagators of the constituents alone. For light quarks these are 
complicated and not so well known functions, however we do know that they have spec- 
tral representations starting from a threshold mass fh 2 and m . Hence the smallest mass 
scale entering the equation comes through this threshold mass. A crude approximation 
of the propagators for p 2 < m 2 and p 2 < m would be, 

S(p) S(p) « -^— 2 (46) 
p — m — m 

where p = p — q. This is qualtitatively reasonable but not quantitatively. In the rest 
frame let us anticipate that there is an average energy p° for quark and p° for antiquark. 
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This can be estimated to be 



m 1 m I A»\ 

Po ~ - - go and q - p « - g (47) 
m + m m + m 



then we have 



S(r>) S(v) ~ - (48) 

Hence naturally in eq. (J48|) the largest length scale given by the exponential fall off of 
the wave function in position space is set by 

5 2 = minim 2 , m)(l - - — ^L— ) (49) 

(m + m) z 

Where \ is the size of the system. For quarks this is an estimate since the propagator is 
not a simple pole. But the existence of the spectral representation for quark propagators 
seems to indicate that it is a good estimate. For deeply bound states the size is entirely 
determined by the threshold mass. The above estimate holds good exactly for long range 
interacting non-relativistic systems such as coulomb interactions. 

This estimate also suggests that when go reaches threshold, 5 vanishes suggesting 
that the bound system attains infinite size. This simple consideration is always valid 
as an estimate whenever there are no massive exchange interactions in the interaction 
kernel. Indeed when there are such interactions the largest length scale between the 
propagators and the interactions (or small mass 5 scale) dominates in dictating the size of 
the system. From the simple minded considerations above we can conclude that massless 
particles cannot be bound as it will necessarily have infinite size. We make an interesting 
observation about chiral symmetric phase of quarks at zero temperature. In this phase 
the quark has vanishing threshold mass hence from our considerations it cannot be bound. 
Chiral symmetric vaccum is automatically a non-confining vaccum. 



6 Discussion 

We have discussed a complete relativistic description of bound states and the BS equation 
is reduced to a set of simpler equations. The form of the equations (fTp^ -(P3 j) is valid for 
a general chiral symmetric interaction kernel. Many of our later algebraic simplifications 
is due to the absence of a mass scale in the interaction kernel however our method can 
handle even if there is a mass scale in the interaction kernel. 

We have concentrated mostly on tightly bound systems primarily because the string 
tension a is large in aQCD. Consequently the tight binding approximation is relevant. 
For a range of low quark masses the lightest meson is the pseudoscalar which extrapolates 
all the way to the Goldstone mode. We verified the PCAC result that the pion mass M n 
is related to the renormalised quark mass m as shown in eq. (|4Uj) . 
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On general grounds we find that there are very few stable mesons. This follows entirely 
from the fact that threshold masses for light quarks rh are much smaller than the string 
tension a. The BS equations can be studied to look for complex eigenvalues and thus 
the unstable mesons. We did find several complex eigenvalues numerically to the set of 
equations ([16 )1 - ([25 )1 . We are as yet unable to systematically analyse them. Primary reason 
being that the method of finding the spectrum is necessary but not sufficient. Another 
important necessary condition we can argue is that of the size of the bound states. For 
stable bound states the argument presented earlier is good enough but for unstable bound 
states this needs to be improved upon. 

Another alternative is to invent a formal series solution as done in WC model. In 
principle this method can be applied here as well but the tedium makes the physics 
non-transparent. Our analytic method of computing the tightly bound spectrum (low 
lying tightly bound states when the coupling is large) was applied to WC model and we 
reproduced the known conventional results [13J. 

Fitting experimental data to this model has shown that the threshold mass for u, d or 
s quarks is about the same because of the string tension a being so large. Consequently it 
is easy to anticipate that there can be large flavour mixing. In our model this is next order 
in ^. Consequently we expect corrections are not necessarily small for light quarks. 

There are several normalisation schemes JI] for the BS wave function such as Cutkosky, 
charge, energy normalisation conditions. One of the primary drawbacks here is that all 
known normalisations for relativistic bound states are not positive norms in the standard 
sense. Consequently they are not of much utitlity. However it has been shown that all 
the known normalisations are equivalent [Tlj . 

We have not dealt with heavy quarks for they fall in a different class altogether. In 
this model string tension decreases at larger energy scales. So the value of string tension 
is indeed much smaller for heavy quarks and thus they fall into the loosely bound regime. 
That is the binding energy is much smaller than the rest mass or the threshold mass. 
This is precisely the non-relativistic limit. If we perform the standard non-relativistic 
approximation to the equations (fT5|) - (j23|) . we do get the standard Schrodinger picture 
[To] in momentum space along with spin-orbit interactions. 

aQCD model has many features of QCD as we explicitly emphasised in our papers 
[HE]. Yet we have shown that for renormalised quark mass(current mass) m > .07\/a we 
disobey a well known inequality of the meson spectrum which is understood theoretically 
and valid experimentally, namely in any flavour sector the pseudoscalar is the lightest 
meson. This follows purely from the positivity properties of QCD in the Euclidean for- 
mulation. Analogous positivity property is not valid in our model. But it is interesting 
to note that it is of no consequence for all light quarks(u, d and s). 

A crude estimate suggests that for heavy quarks in our model i.e., for m > 3^/& we 
recover the pseudoscalar mass inequality. If we consider that a is small for heavy quarks 
we do envisage that charm, bottom, top can also be accomodated. This will be discussed 
in a later publication. But we have to bear in mind that there is a range of quark masses 
which will disobey the light pseudoscalar mass inequality and the quarks that occur in 
nature are not in that regime. 
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We have presented the spectrum calculation explicitly for the case where both flavours 
have the same renormalised mass. We can also do these calculations analytically if they 
are unequal. For the deeply bound states that we have considered the effect is small, 
comparable to i corrections. Finally many of our results can be applied in the context 
of technicolour scenarios as well. 
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